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Abstract 

We initiate the complexity theoretic study of the problem of computing the bits of (real) 
algebraic numbers. This extends the work of Yap on computing the bits of transcendental 
numbers like 7r, in Logspace. 

Our main result is that computing a bit of a fixed real algebraic number is in C = NC X C L 
when the bit position has a verbose (unary) representation and in the counting hierarchy when 
it has a succinct (binary) representation. 

Our tools are drawn from elementary analysis and numerical analysis, and include the 
Newton-Raphson method. The proof of our main result is entirely elementary, preferring to 
use the elementary Liouville's theorem over the much deeper Roth's theorem for algebraic num- 
bers. 

We leave the possibility of proving non-trivial lower bounds for the problem of computing 
the bits of an algebraic number given the bit position in binary, as our main open question. In 
this direction we show very limited progress by proving a lower bound for rationals. 
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1 Introduction 



Algebraic numbers are (real or complex) roots of finite degree polynomials with integer coefficients. 
Needless to say, they are fundamental objects and play in important part in all of Mathematics. 

The equivalence between reals with recurring binary expansions (or expansions in any positive 
integral radix) and rationals is easy to observe. Thus computing the bits of fixed rationals is 
computationally uninteresting. However, the problem becomes interesting if we focus on irrational 
real numbers. Computability of such numbers heralded the birth of of Computer Science in Turing's 
landmark paper [15] where the computability of the digits of irrationals like 7r, e is first addressed. 

Building on the surprising BaileyBorweinPlouffe (BBP) formula [5] for tt, Yap |18j shows that 
certain transcendental numbers such as tt have binary expansions computable in a small complexity 
class like deterministic logarithmic space. Motivated by this result we seek to answer the corre- 
sponding question for algebraic numbers. The answers we get turn out to be unsatisfactory but 
intriguing in many respects. In a nutshell, we are able to show only a very weak upper bound to 
the "succinct" version of the problem and virtually no lower bounds. This gap between best known 
(at least to our knowledge) upper and lower bounds easily beats other old hard-to-crack chestnuts 
such as graph isomorphism and integer factorization. 

1.1 Versions of the Problem 

The problem as stated in [18] asks for the n-th bit of the (infinite) binary sequence of an irrational 
real given n in unary. The succinct version of the problem asks for the n-th bit, given n in binary. 
This version of the problem is naturally much harder than the "verbose" version. We can solve 
the verbose version in C=NC 1 , a subclass of logspace. For the succinct version of the problem 
we are unable to prove a deterministic polynomial or even a non-deterministic polynomial upper 
bound. The best we can do is place it at a finite level in the counting hierarchy (which includes 
the computation of the permanent at its first level) . Even more surprising is that we are unable to 
prove any non-trivial lower bound for any irrational algebraic number. Intriguingly, we can prove 
Parity and in general AC°[p] lower bounds for computing specific rationals. 

1.2 Previous Proof Techniques 

In his article |18| , Yap used a BBP like [I] series to prove a logspace upper bound for the (verbose 
version of) computing the bits of tt. At the core of that argument is the concept of bounded irra- 
tionality measure of tt which intuitively measures how inapproximable tt is, by rationals. Roughly, 
the BBP-like series was used to approximate tt by rationals and then argue, via the bounded irra- 
tionality measure, that, since there aren't too many good approximations to tt, the computed one 
must match the actual expansion to lots of bit positions. 

1.3 Our proof technique 

Further progress was stymied by the extant ignorance of BBP like series for most well-known 
irrationals. Our crucial observation is that approximating an irrational can be accomplished by 
means other than a BBP-like series for instance by using Newton-Raphson. Bounded irrationality 
measure for algebraic numbers follows by a deep theorem of Roth [13]. But we show that we can 
keep our proof elementary by replacing Roth's theorem by Liouville's Theorem [8] which has a 
simple and elementary proof (see e.g. [14]). 
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An upper bound on the succinct version of the problem follows by observing that Newton 
Raphson can be viewed as approximating the algebraic number by a rational which is the ratio of 
two Straight Line Programs or SLP's. Allender et al. [2\ show how to compute bits of a single SLP 
in the Counting Hierarchy. We extend their proof technique to solve the problem of computing a 
bit of the ratio of two SLP's in the Counting Hierarchy. 

1.4 Related Work and Our Results 

Yap [18J showed that the bits of tt are in Logspace. This was the origin of this endeavour and we 
are able to refine his result to the following: 

Theorem 1 Let a be a real number with bounded irrationality measure, which can be expressed as 
a convergent series and further the the m th term (for input m in unary) is in FTC . Then the n th 

ppPP 

bit of a can be computed by a TC circuit for n in unary and in PH for n in binary. 

In particular, the above inclusions hold for tt. 

Somewhat paradoxically we get slightly weaker bounds for algebraic numbers. As our main 
result, we are able to show that (for an explanation of the complexity classes used in the statement 
please see the next section): 

Theorem 2 Let p be a fixed univariate polynomial of degree d, having integer coefficients. Then, 
n th bit of each real root of p can be computed in C=NC 1 nTCLL, if n is given in unary, and in 
ppPP 

PH if n is given in binary. 

Roughly twenty five years ago, Ben-Or, Feig, Kozen and Tiwari [7] studied the problem of finding 
the roots of a univariate polynomial of degree n and m bit coefficients, under the promise that all 
roots are real. Under this assumption they are able to show that approximating the roots to an 
additive error of for an integer \x (which is presumably specified in unary) is in NC. Notice 
that while their result concerns non-constant algebraic numbers it does not involve finding the bits 
of the algebraic numbers only approximating them. Also, since they only achieve unary tolerance 
their claimed upper bound of NC is not better than the C=NC 1 C NC. Also their method does not 
work if the all-real-roots promise is not satisfied. 

1.5 Organization of the paper 

In Section [2] we start with pointers to relevant complexity classes and more importantly known 
results from elementary analysis that we will need in our proofs. In Section [3] we provide upper 
bounds on the complexity of composing bivariate polynomials. This is used in the subsequent 
Section 0] where we view Newton- Raphson as an iterated composition of bivariate polynomials. In 
this section we prove that the method converges "quickly" if its initial point is in a carefully picked 
interval and that we can efficiently identify such intervals. In Section [5] we make use of the tools 
we have put together in the previous sections to prove Theorem [H [2j We also prove a lower bound 
on rationals in this section. Finally in Section [6] we conclude with some open questions. 
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2 Preliminaries 



2.1 Complexity Theoretic Preliminaries 

We start off by introducing straight line programs. An arithmetic circuit is a directed acyclic graph 
with input nodes labeled with the constants 0, 1 or with indeterminates X\, . . . , for some k. 
Internal nodes are labeled with one of the operations +,—,*. A straight-line program is a sequence 
of instructions corresponding to a sequential evaluation of an arithmetic circuit. We will need to 
refer to standard complexity classes like NP, PP, PH and we refer the reader to any standard text 
in complexity such as [3]. We will also use circuit complexity classes like TC°,GapNC 1 and we refer 
the reader to [16] for details. 

One non-standard class we use is TCLL. This is inspired by the class FOLL introduced in [6] 
which is essentially the class of languages accepted by a uniform version of an FAC° circuit iterated 
O(loglogn) many times with an AC°-circuit on top. We obtain a TCLL-circuit by adding a TC°- 
circuit on top of the iterated block of FAC°-circuits. The class of languages accepted by such circuits 
constitutes TCLL. 



2.2 Mathematical Preliminaries 

In order to upper bound the largest magnitude of roots of a polynomial we note the following (e.g. 
see Chapter 6 Section 2 of [T7]): 

Fact 1 (Cauchy) Let p(x) = Yli=o a i xi ^ e a polynomial. Then every root of p(x) is smaller in 
absolute value than: 

M p = l + -, — rmaxflo 1 !^! 

\ a d\ 

We can consider [— M p , M„] as possible solution range which contains all the real roots. 

Fact 2 The Taylor (see [1]) series of a real (complex) function f(x) that is infinitely differentiable 
in a neighborhood of a real (complex) number a is the power series 

ra=0 

where f <yn \a) denotes the n th derivative of f evaluated at the point a. 

We will need to lower bound the minimum distance between the roots of a polynomial or the 
so called root separation. We use the following version of the Davenport-Mahler theorem (see e.g. 
Corollary 29 (i) of Lecture VI from Yap [17] for details of notation and proof): 

Fact 3 (Davenport-Mahler) The separation between the roots of a univariate polynomial p{x) of 
degree d is at least: 

V^s~cjp)\\\p\\ 2 d+l d-^/ 2 

Here, disc(p) is the discriminant of p and ||p||2 is the 2-norm of the coefficients of p. Further 
notice that a lower bound approximation to this bound can be computed in FTC for constant 
polynomials p (since it involves computing lower bound approximations for radicals and powers of 
constants and cosntant determinants). 

We will also need an upper bound on the magnitude of the derivative of a univariate polynomial 
in open interval (a,b). This is given by the so called Markoff's Theorem |llj (see also |12j): 
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Fact 4 Let p(x) be a degree d polynomial satisfying, 

Vxe (a,b),\p{x)\ <M {a>b) 

Then the derivative p'(x) satisfies: 

d 2 M (n M 

y x e(a,b),\p'(x)\<M' ia>b) = -^ 

2.2.1 Removing Rational and Repeated Roots 

The rational roots of a polynomial can be dealt with the aid of the following, easy to see, fact: 

Fact 5 The rational roots of a monic polynomial (i.e. highest degree coefficient is 1) with integer 
coefficients are integers. 

To make p(x) = Ya=o a i x% m onic, just substitute y = adX in a d ^ l p{x) to obtain a monic polynomial 
q(y). Iterating over all integers in the range given by Fact Q] we can find and eliminate the integer 
roots of q(y) and therefore the corresponding rational roots of p{x). 

In general, a polynomial will have repeated roots. If this is the case the repeated root will also 
be a root of the derivative. Thus it suffices to find the gcd g of the given polynomial p and its 
derivative p' , since we can recursively find the roots of p/g and g. As we will be focusing only on 
fixed polynomials the gcd computation and the division will be in FTC . 

We will actually need a more stringent condition on the polynomials - i.e. p does not share a 
root with its derivative p' and even with its double derivative p" . Notice that the above procedure 
does not guarantee this. For example if p(x) = x(x 2 — 1), then p'(x) = 3x 2 — 1 and p"(x) = 6x. 
Thus, p and p" share a root but p and p' don't. 

So we will follow an iterative procedure in which we find the gcd of p, p' to get two polynomials 
Pi = p/(p>p'),P2 = {PiP 1 ) (denoting gcd of p, q by (p,q)). For p\ we find the gcd (pi,p'{) and set 
P3 = Pi/ {Pi->Pi)->P4 = {PiiP'i)- Now we recurse for the 3 polynomials P2,Pz-,P4 (whose product is 
p). Notice that the recursion will bottom out when some gcd becomes a constant. As a result of 
the above discussion we can assume hereafter that p does not share a root with either p' or p" . 

2.2.2 Good Intervals 

Definition 3 Fix an interval I = [a, b] of length \I\ = b — a. We will call the interval I good for 
an integral polynomial p if it contains exactly one root (say a) of p and no root of p' and p" . 

Lemma 4 If p is a polynomial of degree d such that p doesn't share a root with its derivative and 
double derivative then there exist 5 such that all intervals of length less than 5 contain at most 
one root of p and if they contain a root of p then they do not contain any roots of p' or p" . As a 
consequence we can find good intervals I in FTC . 

Proof: Using Fact [3] compute 5 i.e. the minimum distance between roots of p and that of p'p" . (let 
PI1P2 be p',p" with the roots common with p",p"' respectively, removed; then a lower bound on 
the root separation of PP1P2, does the job). Partition the interval containing all the roots (which 
can be inferred from Fact [T]) into sub-intervals of length 5. If such an interval contains a root of p 
then it contains precisely one root of p and no roots of p' or p" . By checking that signs of p are 
opposite at the end points we can identify good intervals. ■ 
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2.3 From approximation to exact computation 



We will crucially use the following theorem (see e.g. Shidlovskii|14j or Yap |17j ) 

Fact 6 (Liouville's Theorem) If x is a real algebraic number of degree d > 1, then there exists a 
constant c = c(x) > such that the following inequality holds for any a £ Z and j3 6 N, a/(3 ^ x: 



a 

x -0 



c 

> ]]d 



The rest of this subsection is an adaptation of the corresponding material in Chee Yap's paper on 
computing ir in L. The primary difference being that we choose to pick the elementary Liouville's 
Theorem for algebraic numbers instead of the advanced arguments required for bounding the irra- 
tionality measure of it. We could throughout replace the use of Liouville's theorem by the much 
stronger and deeper Roth's theorem but prefer not to do so in order to retain the elementary nature 
of the arguments. 

Definition 5 Let x be a real number. Let {x} = x — [x\ be the fractional part of x. Further, let 
{x} n = {2 n x} and x n be the n-th bit after the binary point. 

It is clear that x n = 1 iff {x} n -± > |. For algebraic numbers we can sharpen this: 

Lemma 6 (Adapted from YapfTS\j) Let x be an irrational algebraic number of degree d and let 
c = c(x) be the constant guaranteed by Liouville's theorem. Let e n = c 2 _ ( rf ~ 1 ) n ~ 2 then for n such 
that e n < \ we have: 

• x n = 1 iff {x} n -i G (i + 2e n , 1 - 2e n ). 

• x n = iff {x} n -i £ (2e n , \ - 2e„) . 

Proof: Taking f3 = 2 n in Liouville's theorem we get: \x — 2~ n a\ > i.e., \2 n ~ 1 x — ^1 > 



2 d( n _l) + l 

Consequently, we can find successive approximations {S'mjmgN such that the error terms R r , 
x S m are small enough (as described below). x n is just the n-th bit of S m . 



3 Complexity of Composition 

We investigate the complexity of composing polynomials. This will be useful when we use the 
Newton- Raphson method to approximate roots of polynomials since Newton-Raphson can be viewed 
roughly as an algorithm that iteratively composes polynomials. 

Definition 7 A univariate polynomial with integer coefficients is, an integral polynomial. Any 
integral polynomial when evaluated on a rational value ^, where a, f3 are integers, can be expressed as 
the ratio of two bivariate polynomials in a, (3 called the ratio polynomials of the integral polynomial. 

Definition 8 Let p be an integral polynomial. For a positive integer t, the t-composition of the 
polynomial denoted by p^ is defined inductively as: p^(x) = p(x) and (x) =p it \p{x)). 
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Definition 9 Let f,g be a pair of bivariate polynomials. For a positive integer t define the t- 
bicomposition of(f,g) to be the pair of bivariate polynomials (F^,G^) as follows: 

F®(a,P) = f{a,P),C^{a,P) = g(a,/3), 

and, 

F^(a,ft) = f(F^(a,f3),G^(a,(3)), 
G^(a,(J)=g(FW(a,(3),GW(a,P)), 

The following is a direct consequence of the definitions: 

Proposition 10 The ratio polynomials of the t-composition of an integral polynomial p are exactly 
the t-bicompositions of the ratio polynomials of p. 

Definition 11 For an arithmetic complexity class C containing FTC we call an integral polynomial 
C-computable if its ratio polynomials (viewed as functions that take in the bit representations of its 
two (constant) arguments as inputs and output an integer value) are inC. 

Here we consider upper bounds on the complexity of computing compositions of fixed integral 
polynomials. We first prove that: 

Lemma 12 Let p be a fixed integral polynomial, then given n in unary the l-bicomposition (for 
I = 0(\logn})) of p is computable in GapNC 1 . 

Proof: From Definition [TTJ and Proposition [101 it suffices to prove that the bicomposition of 
the ratio polynomials of p is computable in GapNC 1 . Given a rational number as the bits of 
its numerator and denominator we can first obtain the arithmetic values of its numerator and 
denominator in GapNC 1 . Now we are done modulo the following claim. ■ 

Claim 13 Let f be a fixed bivariate polynomial with integral coefficients and let a, f3 be two integers, 
then there is a constant depth arithmetic circuit that takes a, ft as inputs and outputs f(a,(3). 

Proof: (of Claim) The depth of the circuit is seen to be bounded by 0( [log d~\ ) where d is the degree 
of the polynomial - we just need to find a l f3 J by a tree of height max( [log i], [log j]) + 1 multiply 
it with the coefficient and add up the results by a tree of depth [log (d + 1)] . Since degree is a 
constant we are done. ■ 

We now prove an orthogonal bound on the complexity of compositions. But before that we 
need a small lemma: 

Lemma 14 Suppose G is a layered graph of width 0(n) and depth O(logn) then reachability (from 
a vertex in the first layer to a vertex in the last layer) can be done by an AC-circuit of depth 
0(log logn) 

Proof: It suffices to show that the reachability between two layers which are separated by another 
layer is in AC , which is clear. ■ 

Note 1 In fact, from the proof it is clear that, if G contains O(logn) identical layers then this 
reachability is in the class FOLL defined in f^. 
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Lemma 15 Let p be a fixed integral polynomial, then given n in unary the l-bicomposition (for 
I = 0(\logn})) of p is computable in TCLL. 



Proof: Notice that if a, (3 are some fixed integers, the value of the ^-bicomposition of the ratio 
functions on a, f3 is bounded by an n ^ bit integer. It suffices to compute the value of this 
composition modulo all O(logn) bit primes, since we can do Chinese Remaindering in TC°. Fix 
an O(logn) bit prime q and construct the following bipartite graph H q on vertices S,T (where 
|5| = \T\ = q 2 ) and both S, T consist of pairs ab, a,b 6 {0, . . . ,q — 1}. (ab,a'b ! ) is an edge iff 
f(a, b) = a' mod q and g(a, b) = b' mod q where /, g are the ratio functions of p. Further, let G q be 
the layered graph obtained by taking I layers of H q . Then it is clear that reachability in G q from 
the first to the last layer is exactly equivalent to the values of the Z-compositions modulo q. We are 
done with the aid of Lemma [14] and [9] . ■ 

The following is a consequence of the definitions and of [9]. 

Lemma 16 If p is an integral polynomial which is C-computable (FTC°C C), then on input m,n 
in unary, where m > n, we can obtain the n-th bit of some number that differs from p(^), by at 
most 2~ m in C. 

Now we describe the binary analog of the above lemma. 

Note 2 In the remaining part of this section we denote polynomially bounded integers by lower 
case letters e.g. n,t. We denote those with polynomial number of bits by uppercase letters e.g. 
N,T. Finally we denote those with exponentially many bits by calligraphic letters e.g. ftf,T>. This 
notation does not apply to rationals like u,o~. 

Lemma 17 Let Af and T> be the outputs of two SLP's. Computing the N th (where N is input in 



binary) bit of an approximation (accurate up to an additive error of 2 ( JV+1 ) ) of £j is in PH 



PP 



PP 



Proof: We will compute an under approximation of ^ with error less than 2~( N+1 \ 

Let u = 1 — V2~ T where T > 2 is an integer such that 2 T ~ X <T>< 2 T . Hence |n| < i. 
Notice that the higher order bit of T can be found by using PosSLP : we just need to find an 
integer t such that 2 2 < T> < 2 2 and both these questions are PosSLP questions. Having found 
Ti a lower bound of T correct up to the higher order i bits of T, i.e. 2 Ti < T> < 2 Ti+21 , we check 
if 2 Ti+21 < T> and update Tj_i to T, + 2* _1 iff the inequality holds (and Tj_i = Tj otherwise). 
Thus by asking a polynomial number of PosSLP queries, we can determine T, so each bit of T is 

ppPP 
in PH KK . 

Now consider the series 

V- 1 = 2- T (l - u)- 1 = 2- T (l +u + u 2 + ...) 
Set V = 2- T {\ + u + u 2 + ...u N+1 ), then 

£>-!_£>'< 2 - T 2 ' 1 < 2-( W+1 ) 

I>N+1 
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Now we need to compute N th bit of 

K rN+l N+1 
7=0 7=0 

We need to compute the M = N + (N + 2)T th bit of: y = J2i=o~N{2 T - V) I 2 i ~ N+l - 1 ^ '. Since 
each term in summation is large and there are exponentially many terms in summation, so we will 
do computation modulo small primes. Let yj denote the I th term of summation. 

2 

Let A4 n be the product of all odd primes less than 2™ . For such primes P let Hp tU denote 
inverse of ^p 1 mod P. Any integer < Yj < Ai n can be represented uniquely as a list (Yj p), 

2 

where P runs over the odd primes bounded by 2 n and Yi p = 3^7 mod P. 

Define the family of approximation functions app n (y) to be J2p 5^7 Yi,pHp,n&p,n where ap :Tl is 

1 ^ 

the result of truncating the binary expansion of -p after 2™ bits. Notice that for sufficiently large 

3 

n, and y < A4 n , app n {y) is within 2 -2 " < 2~( Ar+1 ) of y /A4 n as in the proof of Theorem 4.2 of 
[2]. Continuing to emulate that proof further and using the Maciel-Therien (see [TO]) circuit for 

ppPP 

iterated addition, we get the same bound as for PosSLP in 2 viz. PH rr . Notice that we have 
a double summation instead of a single one in [2j, yet it can be written out as a large summation 
and thus does not increase the depth of the circuit. ■ 



4 Establishing Quadratic Convergence 

We use the famous Newton-Raphson method to approximate Algebraic Numbers. The treatment is 
tailored with our particular application in mind. There are some features in the proof (for instance 
a careful use of Markoff 's result on lower bounding the derivative of a polynomial) which led us to 
prove the correctness and rate of convergence of the method from scratch rather than import it as 
a black-box. 

Definition 18 (Newton-Raphson) Given an integral polynomial p and a starting point xq, recur- 
sively define: 

_ _ P(xi) 

P{Xi) 

whenever Xi is defined andp'(xi) is non-zero. 

Recall good intervals from Definition [3l 

Definition 19 Given a good interval I for an integral polynomial p, let ei denote the error in the 
v iteration of Newton-Raphson | when starting with xq E I. Notice that is 

defined only when X{ is. 

Definition 20 We say that Newton-Raphson converges quadratically (with parameter M) for an 
integral polynomial p whenever M is a non-negative real such that for any interval I which is good 
for p and of length at most min(^g2",|), it is the case that the errors at consecutive iterations 
(whenever both are defined) satisfy e^+i < Mef. 

The following Lemma shows that not only are the errors at all iterations defined under the 
assumptions of Lemma [20] but also, that, Newton-Raphson converges "quickly". 
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Lemma 21 If Newton- Raphs on converges quadratically (with parameter M) for an integral poly- 
nomial p, then for every i > 0, the i th iterand, Xi, is at distance at most mi\i{^p, 2~ 2 ' /2 ) from the 
unique root of p in any good interval I of length \I\ < min(|, tot)- In particular, Xi £ I for every 
i>0. 

Proof: We proceed by induction on the number of iterations. For the base case, notice that xq is 
at distance at most eo < min(|, j^i) from the root. 
Now assume that e« < m.in(j^ 7 ,2~ 2 '' /2 ). Then, 



< Me 2 



1.5 

i 



= 

2 1 

< 2- 1 2~ L5x2l/2 

2 _ 2 (i+l)/2 



Since e^+i < \e\- 5 and ef 5 < (i j; 1 ) /2 < 1 for i > 0, therefore, e^+i < e{ < where the second 
inequality follows from the inductive assumption. This completes the proof of the inductive step. 



Lemma 22 For any integral polynomial p and any good interval I thereof, there exists a subinterval 
I' C I such that Newton-Raphson converges quadratically in I'. 

Proof: Let the unique root of the integral polynomial p, contained in /, be a. Thus, p(a) = 0. By 
Taylor's series 

p(a) = = p(xi) + (a - Xi)p'(xi) + -(a - Xi) 2 p"(^) 
where £j is between Xi and a. Rearranging, and using the equation for Xi + \ 



a - Xi+i = (a - Xi) + — 



p'(xi) 2p'(xi) 



(a - Xif 



On the other hand, by Definition [19] the error in the i + 1 th iteration of Newton-Raphson (whenever 
defined) is: 



e i+ i = \x i+ i - a\ 



2 p'( Xi ) 

Since p' does not have a root in that interval and p" is finite (because p is a polynomial), the right 
hand side is well-defined. 

In the good interval /, p' is monotonic and hence the minimum (and maximum) value of p' 
is attained at the end-points of /. Now we can upper bound the absolute value \p"\ using upper 
bound for p' and Markoff's result (FactHJ). Let this value be denoted by p± and the minimum value 
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of \p'\ by p2 7^ (p2 = would contradict the assumption that I does not contain a root of p'). 
Now, set M to be ^. 

Partition I into sub-intervals of length and let /' be the unique subinterval containing a 
root of p : i.e. the unique sub-interval such that p takes oppositely signed values at its end points. 
It is easy to see that Newton- Raphson converges quadratically (with parameter M) in I'. ■ 

5 Putting it all together 

We now complete the proofs of Theorem [2] and Theorem [TJ 

Proof:(of Theorem[2D From Lemma[2]we can compute a good interval /. Then using Lemma [221 we 
can find a subinterval of / such that Newton-Raphson will converge quadratically in this interval. 

Since Newton-Raphson converges quadratically, in order to obtain an inverse exponential error 
in terms of n, (By Lemma [2T]) we need 0([~logn~|) iterations. Now by Lemma [121 and [T5l along 
with Lemma [161 we get that 0([~logn~|) compositions of Newton-Raphson (taking as initial point, 
the middle point of the interval /' obtained from Lemma [22]) can be computed in C=NC 1 nTCLL. 
Finally Lemma [6] ensures that we have computed the correct bit value. The argument for the binary 
case is analogous and uses Lemma [T7] instead of Lemmas [121 021 and [161 ■ 

Proof: (of Theorem [I]) Let a be a constant have series of the form a = Yl'kLo = S n + R n where 
we have split the series into a finite sum S n = Efc=o*fc and a remainder series R n = EfcLn+i^fc- 
Each term of series can be written as a rational number of the form if. = /3 -fcc ^y where /3 is a 
real number p(k), q(k) are fixed polynomial with integer coefficients and c > 1 is an integer. 

This series consist of summation of iterated multiplication, division and addition which can be 
computed by TC° circuit. Since a has bounded irrationality measure so its n th bit can be computed 
using Lemma [6l ■ 

Using the BBP-like series for it [I] and its bounded irrationality measure, we get: 

„ ppPP 

Corollary 23 Computing n bit ofn is in TC and PH , given n in unary and binary respec- 

tively. 

5.1 Lower Bound 

Finally we show the Mod p (for any odd prime p) hardness of the bits of a rational. We still don't 
have the proof of any kind of hardness of an irrational algebraic number. 

Lemma 24 For given a odd prime p and an integer X (having binary expansion 6 n _i . . . bo) then 
there exist an integer N , whose bits are constructible by Dlogtime uniform projections, and a fixed 
rational number Q such that N th digit in binary expansion of Q is iff b{ = (0 mod p) . 

Proof: For a given odd prime p we can find a integer t, < t < p such that 2* = l(modp). Such 
a t exists because the multiplicative group of integers modulo p is finite. Consider the number 
N = ^2™=q 6i(2*)* . Then, N = Y17=o &i(modp), because 2* = l(modp). Now, consider the sum: 

Q = Etv>o = g |Swffef The Nth di g it; of Q is iff Ei &i = mod P- ■ 
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6 Conclusion 



We take the first step in the complexity of Algebraic Numbers. Many questions remain. We focus on 
fixed algebraic numbers - in general we could consider algebraic numbers defined by polynomials of 
varying degrees/coefficients. We have ignored complex algebraic numbers - they could present new 
challenges. Most importantly, our study is, at best, initial because of the enormous gap between 
lower bounds (virtually non-existent) and the upper bounds. Narrowing this gap is one of our 
future objectives. 
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